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Abstract 

A phenomenological two-fluid model of the (time-reversible) spectrally-truncated 
3D Euler equation is proposed. The thermalized small scales are first shown to 
be quasi-normal. The effective viscosity and thermal diffusion are then determined, 
using EDQNM closure and Monte-Carlo numerical computations. Finally, the model 
is validated by comparing its dynamics with that of the original truncated Euler 
equation. 



1 Introduction 

It is well-known that the (inviscid and conservative) truncated Euler equation 
admits absolute equilibrium solutions with Gaussian statistics, equipartition of 
kinetic energy among all Fourier modes and thus an energy spectrum E(k) ~ 
k 2 PP . Recently, Cichowlas et al. [2f3] observed that the Euler equation, with a 
very large (several hundreds) spectral truncation wavenumber fc max , has long- 
lasting transients which behave just as those of high-Reynolds number viscous 
flow; in particular they found an approximately k~ 5 ^ 3 inertial range followed 
by a dissipative range. How is such a behavior possible? It was found that the 
highest-/c modes thermalize at first, displaying a k 2 spectrum. Progressively 
the thermalized region extends to lower and lower wavenumbers, eventually 
covering the whole range of available modes. At intermediate times, when the 
thermalized regime only extends over the highest wavenumbers, it acts as a 
thermostat that pumps out the energy of larger-scale modes. Note that similar 
£r 5 / 3 / k 2 spectra have already been discussed in the wave turbulence literature 
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(e.g., [4]) and were more recently obtained, within a simple differential closure, 
in connection with the Leith model of hydrodynamic turbulence [3]. 

The purpose of the present work is to build a quantitative two-fluid model 
for the relaxation of the 3D Euler equation. In section |2j after a brief recall 
of basic definitions, the statistics of the thermalized small scales are stud- 
ied during relaxation. They are shown to be quasi-normal. Our new two-fluid 
model, involving both an effective viscosity and a thermal diffusion, is intro- 
duced in section |3j The effective diffusion laws are then determined, using an 
EDQNM closure prediction and direct Monte-Carlo computations. The model 
is then validated by comparing its predictions with the behavior of the original 
truncated Euler equation. Finally section [4] is our conclusion. 



2 Relaxation dynamics of truncated Euler equations 



2.1 Basic definitions 

The truncated Euler equations ([I]) are classically obtained pQ by performing a 
Galerkin truncation (v(k) = for sup a \k a \ > k max ) on the Fourier transform 
v(x, t) = J2 v(k, t)e ik ' x of a spatially periodic velocity field obeying the (unit 
density) three-dimensional incompressible Euler equations, d{v + (v ■ V)v = 
— Vp , V ■ v = 0. This procedure yields the following finite system of ordinary 
differentials equations for the complex variables v(k) (k is a 3 D vector of 
relative integers (^1,^2,^3) satisfying 

d t v a (k,t) = --P Q/37 (k)^{) / 3(p,t){) 7 (k- p,t) (1) 
1 p 

where V a ^ = kpP ai + k^P a/3 with P al3 = 5 a p — k a kp/k 2 and the convolution 
in (fl]) is truncated to sup a | SUp a \p a \ < k 

max 

and sup a \ k Q — p a \ < 

max- 

This time-reversible system exactly conserves the kinetic energy E = J2k E(k, t), 
where the energy spectrum E(k, t) is defined by averaging v(k', t) on spherical 
shells of width Ak = 1, 

E(k,t) = l E Iv(k',i)| 2 . (2) 

k-Ak/2<\k'\<k+Ak/2 
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2.2 Small Scales Statistics 



Perhaps the most striking result of Cichowlas et al. [2] was the spontaneous 
generation of a (time dependent) minimum of the spectrum E(k, t) at wavenum- 
ber kth(t) where the scaling law E(k,t) = c(t)k 2 starts. Thus, the energy dissi- 
pated from large scales into the time dependent statistical equilibrium is given 
by 

E th (t)= J2 ^(M)- (3) 

k th (t)<k 

In this section we use the so-called Taylor-Green [6] initial condition to ([T|: the 
single-mode Fourier transform of u TG = sin x cosy cos z, v TG = —u (y,—x,z), 
w TG = 0. 

In order to separate the dynamics of large-scale (k < k t h) and the statistics of 
small-scales (k > kth) we define the low and high-pass filtered fields 

/<(r) = £F(k)/ k e lk - r (4) 
/>(k) = l-/<(r) (5) 

where /(r) is an arbitrary field and ff. its Fourier transform; we have chosen 
F(k = § (l + tanh [ lk ^ k ] ) , with A A; = 1 /2. 

This filter allows us to define the large-scale velocity field v < and the spatially 
dependent thermalized energy (or heat) associated to quasi-equilibrium. Using 
the trace of the Reynold's tensor [7j , Rij = ^(vf'vf) < , we define the local heat 

as 

QW = ^[(v>) 2 ] < (r). (6) 

By construction of the filters, (|4]|5]) the heat spatial average is equal to the 
dissipated energy ^ < Q(r) >= £th- Fig. [T^i shows a 2D cut of the heat Q 
on the surface z — \ , where a cold zone is seen to be present at the center of 
the impermeable box (x = [0, 7r], y = [0, 7r], z = [0, 7r]). An isosurface of the 
hottest zones is displayed on Fig. [T|d. Is is apparent on both figures that Q(r) 
is not spatially homogeneous. 



2.3 Heat diffusion 



The simplest quantities to study in order to quantify the evolution of Q, 
are the spatial average Q(t) = (Q(r,t)) and the root mean square variation 
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Fig. 1. Cut at z — \ of Q (a) and the isosurface Q(r) — 0.8Q m ax — 0.42 (b). 



= J ((Q 2 — (Q) 2 ))- These quantities are shown in figure 2 1 where that the 
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Fig. 2. Plots of Q(t) (a) and AQ(t)/Q(t) (b); solid lines are the results of the 
two-fluid model (see section [3]) . 

mean heat is seen to increases in time, due to the energy coming from the 
large eddies, as was shown precedently in |2]. The relative fluctuation AQ/Q 
is seen to decrease from 0.9 to 0.2. 

The next natural question is related to the statistical distribution of the small 
eddies v > : are they approximately Gaussian, like an absolute equilibrium? A 
histogram of v> is shown in figure [3]. As the heat is not homogeneous, we also 
computed the histogram of the normalized field v> — v x I which seems to 
better obey Gaussian statistics as can be seen on figure [3] and comparing the 
firsts normalized cumulant s n = -y= (c n is the cumulant of order n) in the 

table. 



Fig. 3. Histogram of v> and v% an d normalized cumulant S4 and sq (odd cumulants 
vanish because of symmetries). 

3 Two-fluid Model 



We now introduce our phenomenological two-fluid model of the truncated 
Euler equation. One of the fluids describes the large scale velocity field and 
the other represents the thermalized high-wavenumber modes described by 
a temperature field T = Q/c (c is the specific heat, explicitly given by 
c = 8£^ ax ). This model is somewhat analogous to Landau's standard two- 
fluid model of liquid Helium at finite temperature T where there is a natural 
cut-off wavenumber for thermal excitations: the classical-quantum crossover 
wavenumber & max given by hk max cs = ksT (eg is the sound velocity and fc# 
Boltzmann's constant). In Landau's model /c max is temperature dependent and 
the specific heat c is proportional to T 3 . In constrast, k max and the specific 
heat are constant in our model that reads: 



d t vf + vfdjvf = -dip + dp'^ (7) 

div< = (8) 
1 

2c 



d t T + vfdjT = VT + - (djV< + divf) a\ 3 (9) 



where 



cr'ij = F Wes(k) (ihvf + ikjvf)] (10) 
VT = F- 1 [-k 2 D cS (k)F[T]] (11) 

and denotes the Fourier transform, a'^ is a generalized form of the stan- 
dard viscous strain tensor [8] . The precise form of the anomalous diffusion 



terms i/ e g and D c g will be determined below, in sections 3.1 and 3.2 



The advection terms in equation ([7]) are readily obtained from the Reynolds 
equations for the filtered velocity by remarking that the diagonal part of the 
Reynolds stress can, because of incompressibility, be absorbed in the pressure. 
Equation ( fIo| represents a simple model of the traceless part of the Reynolds 
tensor [7j. In the same vein, the advection terms in equation (|9j) are readily 
obtained together with higher-order moments (see equation (1) of reference 
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[9]). The dissipation and source terms in ^ are thus simple models of the 
higher-order moments. It is easy to show that in the present model (|v 2 + cT) 
is conserved, corresponding to the energy conservation in the truncated Euler 
equation. 

As the fluctuations AQ/Q are small (see above) we will furthermore assume 
that u eS and D eS only depend on < Q >= E th . Thus the evolution of the 
filtered velocity v< is independent of the fluctuations AQ. As [-Eth] = L 2 T~ 2 , 
simple dimensional analysis yields the following form for the function v c & and 
D eS : 

y/E^ f f k k \ >/^h lT , ( k k \ . . 

" eS = Jc f [k~'k~) ;jDcff = l J (12) 

"Triax \ "'max "rnax / "'max \ "'max "'max / 

where k = 2tt/L p the smallest nonzero wavenumber (L p is the periodicity 
length, 2tt in the present simulations). 



3.1 ED QNM Determination of viscosity 



An analytical determination of function v eS is possible using the eddy-damped 
quasi-Markovian theory (EDQNM) [TO] . It is known that this model well re- 
produces the dynamics of truncated Euler Equation, including the k~ 5 ^ 3 and 
k 2 scalings and the relaxation to equilibrium [TT] . 



The EDQNM closure furnishes an integro-differential equation for the spec- 
trum E(k, t): 

dE(k,t) 

— — — = T NL {k,t) (13) 
where the nonlinear transfer Tjvx is modeled as 



T NL {k,t)= J J Q kpq {xy + z 3 )[k 2 pE{p,t)E{q,t) 



A 

3 , j.Mn/i ^dpdq 



-p d E(q,t)E(k,t)}^. (14) 
pq 



In (14) A is a strip in p, q space such that the three wavevectors k, p, q form 
a triangle, x, y, z, are the cosine of the angles opposite to k, p, q. Qk pg is a 
characteristic time defined as 

n 1 -exp(-(Vk + V P + Vq)t) n _v 

^kpq — \±0) 
Vk + Vp + Vq 
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and the eddy damped 7] is denned as 



Vk 



\ 



s 2 E(s,t) ds. 



(16) 



Classically A = 0.36 and the truncation is imposed omitting all interactions 



involving waves numbers larger than /c max in (14). 



A simple and important stationary solution of (13) is the absolute equilibrium 



with equipartition of the kinetic energy and corresponding spectrum E(k) ~ 
k 2 . 

To compute the EDQNM effective viscosity u e s we consider an absolute equi- 
librium with a small perturbation added in the mode of wavenumber k pert and 
study the relaxation to equilibrium. The corresponding ansatz is E(p, t) = 
W^p 2 + r y{t)5{p — A; pcrt ) and we suppose £7 t h 3> 7, so that the total energy is 
almost constant and equal to E^. 



Using the long time limit of (15) and expanding the EDQNM transfer (14) 



to first order in 7 yields for the delta containing part, after a lengthy but 
straightforward computation: 



T NL (k,t) = -y(t)5(k - k peTt )k 



2 v E th y 30 



k 



(17) 



where I is given by the explicit integral 



I{x) = \fx 



(p 2 -l)(l-g 2 )(g 2 +p 2 (l + 2g 2 ))) 
' (p 2 -g 2 )(2§ +((p- g )f +(p + g )f)) 



dqdp 



Using (13), (17) and the basic definition of the two-fluid model ([7 11), we 
obtain 



v e s{k) 



^Vsof k 



ki] 



k 9\ 

The function f(x = ^r^, 0) in (12) is thus given by 
7 30 



f(?,0) 



2A 



I x). 



In the limit x — > 0, it is simple to show that / has a finite value /(0, 0) 
Thus the EDQNM prediction in the small k/k max limit is 



(19) 



15A' 



^cff 



th 



15A' 



(20) 
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with 



'15A 



5.021 for the classic value of A = 0.36. This asymptotic value can 



also be obtained from the EDQNM eddy viscosity expression calculated by 
Lesieur and Schertzer [12] using an energy spectrum E(k) ~ k 2 . 



3.2 Monte- Carlo determination of viscosity and thermal diffusion 



In order to numerically determine the effective viscosity v e s{k) of the two- fluid 
model, we use a general-periodic code to study the relaxation of an absolute 
equilibrium perturbed by adding a stationary solution of the Euler equation. 
We thus consider the initial condition 



u = cos kx sin ky + u cq 
v = — sin kx cos ky + v eq 



w = w. 



cq 



(21) 
(22) 
(23) 



where the (solenoidal and Gaussian) absolute equilibrium velocity field satis- 



fies « + ^ + = 2E- 



th- 



The resulting amplitude of the rotation in (21 23) is found, after a short 



transient, to decay exponentially in time. The function v e s(k) is then obtained 




b) 
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Fig. 4. Effective viscosity u c g (a) and thermal diffusivity D e g (b) determined by 
Monte Carlo computations performed at different values of E t h and /c max (see text). 

by finding the halving time Tfc, for which u a (k, to + T&) = u a (k, to)/2, with to 
chosen larger than the short transient time. The effective dissipation thus 
reads 



v eS (k) = log2/(fc 2 r fc ). 



(24) 



The values of v e s(k)k ma/X / ^E th are shown in figure (|4p,) for different values of 
E th , k, fc max . A very good agreement with the EDQNM prediction is observed. 
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Note that there is not dependence in the dimensionless parameter /co//c max (see 



eq. Ql2p). 

An exponential fit of all data in figure |4^i gives 



z/ eff = 5.0723^e- 397fc / fc " 



k'y, 



Note that the limit k/k n 



(25) 



is consistent with the EDQNM prediction 



(20) 



Another simple numerical experiment can be used to characterize the thermal 
diffusion: the relaxation of a spatially-modulated pseudo-equilibrium defined 
by 



(u 2 + v 2 + w 2 ) = 2E th + 2e cos (fee) 
with e < E th . 



(26) 



An x-dependent temperature can be recovered by averaging u 2 + v 2 + w 2 over 
y and z. Numerical integration of the truncated Euler equation with the initial 



condition (26) produces an amplitude e that decays exponentially, as in the 



case studied for the determination of effective viscosity The thermal diffusivity 



D e ?( is determined in the same way as in eq. (24) and the corresponding data 
is shown in figure [4b. A power-law fit gives 



D 



,,r = 0.7723^(A;/A; n 



-0.74 



(27) 



The negative exponent in (27) is characteristic of hypodiffusive processes. 




Fig. 5. Effective Prandtl number P e R = v e g/D c g. 



We can define an effective Prandtl number as the ratio P c s(k) = v e s(k) / D e g(k). 
The Prandtl number is plotted in figure [5j where the solid blue line is obtained 
using the EDQNM prediction (20) and the fit (27) and the dashed red line is 



obtained using the fits (25) and (27). Note that the Prandtl vanishes in the 



the small k/k max limit and verifies P e ff < 1 for all wavenumbers. 
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3. 3 Validation of the Model 



In this section, numerical integration of the the two-fluid model equations (J^- 
11) are performed using a pseudo-spectral code. Time marching is done using 
second-order leapfrog finite difference scheme and even and odd time-steps are 
periodically re-coupled by fourth-order Runge-Kutta. The effective viscosity 
and diffusivity are updated at each time step by resetting E t ^ =< Q >. The 
obtained data is compared with that directly produced from the truncated 
Euler equation. 




Fig. 6. a)Time decay of rotation (21 22]) (upper curve) and temperature modulation 
(26) (bottom curve). Solid line: truncated Euler equations and dashed line: two-fluid 
model, b) Time-evolution of energy spectra, truncated Euler equation: solid lines and 
two-fluid model: dashed lines. 



The time-evolutions resulting from initial data (21 22) (in red) and p6| ) (in 
blue), both normalized to one and with the same value of E t ^ is displayed 
on figure |6^ . Good agreement with the two-fluid model is obtained in both 
cases and the faster relaxation of the temperature modulation is related to 
the smallness of P e g < 1- 

We now compare, the evolution of non-trivial spectra of the truncated Euler 
equation and the two-fluid model. The truncated Euler equation is in- 
tegrated using the Taylor-Green initial data. At t ~ 8, when a clear scales 
separation is present, the large-scale fields v < (see eqjij and the heat Q (eq|6]) 
are computed and used as initial data for the two-fluid model ((7 11). The 



subsequent evolution of the two-fluid model is then compared with that of the 
truncated Euler equation. 

Both spectra, plotted in figure [6b, are in good agreement. The straights lines 
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represents the thermalized zone E(k,t) = c(t)k 2 in the the spectrum of the 
truncated Euler equation, where c(t) is determined by the condition (Q(t)) = 

£ c{t)k\ 

fc>fc t h 

The value of Q(t) and AQ/Q are plotted in figure [2] (solids lines); the evolution 
of the fluctuation of the temperature are well reproduced too by the two-fluid 
model. 



4 Conclusion 



The thermalized small scales were found to follow a quasi-normal distribu- 
tion. The effective viscosity was determined, using both EDQNM and Monte 
Carlo. (Hypo)diffusion of heat was obtained and the effective Prandtl num- 
ber found to vanish at small fc/fc max . The two-fluid model was found to be in 
good quantitative agreement with the original truncated Euler equations. We 
acknowledge useful scientific discussions with A. Pouquet and the support of 
ECOS and CONICYT. 
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